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The  objectives  of  this  research  were  to  develop  and  refine  a  shape  memory  alloy  (SMA) 
modeling  approach,  the  martensite  twin  reorientation  (MTR)  snbrontine,  to  predict  the 
transient  response  of  a  spaceborne  membrane  optic  SMA  deployment  actnator  spine  system. 
In  concert  with  a  commercial  finite  element  solver,  this  application  snpports  mnlti- 
dimensional,  temperatnre-displacement  transient  predictions  of  the  shape  memory  effect 
(SME)  exhibited  by  SMAs  throngh  implementation  of  a  phenomenological  constitntive  law. 
The  scope  of  this  stndy  was  to  model  the  relation  of  inpnt  power  magnitnde  and  waveform  to 
stress  fields,  reaction  forces,  and  thermal  fields  for  the  fignre  acqnisition  of  a  gossamer 
reflector.  Of  particnlar  interest  is  the  stress  and  thermal  field  history  of  the  polymeric 
membrane  concentrator  throngh  the  deployment  to  prevent  mechanical  and  thermal  failnre 
as  limited  experimental  or  modeling  analysis  resnlts  exists  for  SMA  deployment  schemes. 
Thermal  anthority  over  the  concentrator  was  fonnd  locally  limited  to  the  spine  interface. 
Upon  apertnre  engagement,  the  integrated  model’s  first  excnrsion  cycle  freqnency  excited  by 
all  three  loading  regimes  was  non-nniqne,  regardless  of  the  preceding  temporal  character  of 
the  spine’s  response.  This  tendency  snggests  a  conpling  between  the  strnctnral  dynamic 
designs  for  deployment  and  in-service,  deployed  architectnres.  Aspects  of  the  deployment 
seqnence  warranting  fnrther  stndy  and  issnes  for  tackle  to  farther  develop  the  MTR 
snbrontine  are  identified.  Technology  addressed  throngh  this  thesis  research  is  intended  to 
foster  and  matnre  snccessive  large,  lannch-packaged  space  vehicle  programs. 
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Nomenclature 

B2  phase  finish  temperature 
B2  phase  start  temperature 
concentrator  specific  heat 
spine  specific  heat 
elastic  modulus 
B2  phase  modulus  of  elasticity 
concentrator  modulus  of  elasticity 

modulus  of  elasticity  based  upon  element  number  i  and  loop  number  j 

R  phase  modulus  of  elasticity 

spine  modulus  of  elasticity 

focal  length  to  diameter  ratio 

concentrator  focal  length  to  diameter  ratio 

spine  focal  length  to  diameter  ratio 

force 

nodal  B2  phase  geometry  vector  based  upon  node  number  i 
nodal  geometry  vector  based  upon  node  number  i  and  loop  number  j 
nodal  R  phase  geometry  vector  based  upon  node  number  i 
concentrator  thickness 
spine  thickness 
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I.  Introduction 

SINCE  discovery,  the  shape  memory  effect  (SME)  or  martensite  twin  reorientation  (MTR)  has  been  identified  as 
a  potential  agency  to  drive  active  material  systems.  Shape  memory  alloys  (SMAs)  are  particularly  attractive  for 
space  applications  due  to  their  ability  to  recover  large  deformations  by  generating  high  specific  force,  with 
reasonable  energy  consumption  and  intrinsic  precision  movement,  an  artifact  of  the  atomistic  nature  of  the  austenite- 
martensite  crystalline  phase  transformation.  Although  considerable  advancements  have  been  made  in  the  evolution 
of  tunable,  reconfigurable,  and  actuation  material  systems,  the  crux  for  progress  of  this  technology  is  our  ability  to 
model  the  thermo-mechanics  of  the  SME.* 


A.  Background 

SMA  components  first  established  flight  heritage  on  space  vehicles  as  alternatives  to  pyrotechnic  devices. 
Applications  include  the  back-up  boom  release  on-board  the  ISEE-B  Spacecraft  (1977),  a  solar  array  bearing  pin 
off-load  and  unlatching  mechanism  for  the  Hubble  Space  Telescope  solar  panels  (1990),  and  a  solar  panel  release 
on-board  the  Clementine  Spacecraft  (1994).^  Just  as  these  vehicles  incorporated  the  SME,  deployment  and  in- 
service  challenges  of  gossamer  spacecraft  will  further  showcase  the  leverage  SMA  technology  extends  to  structural 
reconfiguration  and  tuning. 

Gossamer  technology  is  anticipated  to  enable  a  paradigm  revision  in  the  approach  to  space  system  design 
philosophy  through  revolutionary  spacecraft  component  technologies  and  driven  by  mission  requirements 
demanding  order  of  magnitude  increases  in  aperture  size  for  higher  spatial  and  temporal  resolution.  Challenges  of 
this  emerging  revision  are  deployment  techniques,  membrane  and  rigidizable  materials,  surface  accuracy  for 
precision  reflectors,  and  analytical  tool  development.  Pioneering  gossamer  spacecraft  include  the  Echo  balloon 
series  (1950s  and  1960s),  Goodyear  antennas  (1960s),  Contraves  antennas  and  sunshades  (1970s  and  1980s), 
L’Garde  inflatable  decoys  (1970s  and  1980s),  and  L’Garde  Inflatable  Antenna  Experiment  (1996).^  Technology 
addressed  through  this  thesis  research  is  intended  to  foster  and  mature  successive  large,  launch-packaged  space 
vehicle  programs. 

B.  Literature  Review 

The  earliest  documentation  of  the  SME  is  attributed  to  Chang  and  Read"*  in  1951.  Buehler  and  Wiley,^  in  1962, 
observed  this  phenomenon  in  a  popular  SMA  while  conducting  materials  studies  at  the  Naval  Ordinance  Laboratory. 
This  particular  equiatomic  nickel-titanium  alloy  was  identified  as  Nitinol  to  incorporate  the  namesake  of  both  the 
elemental  composition  and  laboratory  of  discovery.  Since  Buehler  and  Wiley’s^  find,  considerable  research  effort 
has  been  invested  to  understand  the  micro-structural  characteristics  involved  and  to  mathematically  describe  the 
response  of  Nitinol,  as  well  as  other  metallic  alloys  found  to  exhibit  similar  behavior,  such  as  the  Cu-Al-Ni  and  Cu- 
Zn-Al-Mn  systems.  The  considerable  attendant  recoverable  inelastic  strains  (8-9  %)  and  stresses  (500-900  MPa  for 
Nitinol)  of  the  SME  have  perennially  sustained  interest  in  this  materials  technology  amongst  the  aerospace 
community.^’  ’ 

Several  constitutive  models  have  been  proposed  to  mathematically  describe  the  macroscopic  mechanical 
behavior  of  SMAs.  Most  notably  among  those  considered  a  phenomenological  approach  include  constitutive  laws 
developed  by  Tanaka  and  Nagaki,*  Liang  and  Rogers,^  Ivshin  and  Pence,'®  Brinson,"  Boyd  and  Lagoudas,'^  and 
Abeyaratne  et  al.'^  These  constitutive  laws  generally  incorporate  both  a  kinetic  law  governing  the  phase 
transformation  and  a  mechanical  law  governing  stress-strain  behavior.  Often  the  kinetic  laws  distil  to  a  single  state 
variable,  the  martensite  fraction,  as  a  function  of  temperature  and  stress.  Models  employing  these  laws  yield 
relatively  accurate  uni-axial  solutions,  as  they  are  founded  on  ad  hoc  descriptions  sought  to  agree  with  empirical 
results.  The  most  significant  discrepancy  between  these  models  is  the  imbedded  kinetic  law  as  the  detwinned  strain 
often  dominates  the  conventional  elastic  deformation  mechanism.’"' 

Models  based  on  a  micro-mechanics  approach,  founded  on  thermodynamic  principles,  include  those  developed 
by  Sun  and  Hwang,’^  Patoor  et  al.,’®  Goo  and  Lexcellent,’^  Huang  and  Brinson,’*  Vivet  and  Lexcellent,’®  and  Lu  and 
Weng.^®  Theses  micro-mechanical  formulated  constitutive  laws  also  yield  relatively  accurate  uni-axial  solutions, 
but  are  more  computationally  demanding  due  to  increased  complexity.  A  mesomechanical  model  has  also  been 
proposed  by  Vokoun  and  Kafka^’  based  on  the  general  concept  deformation  of  two-phase  materials  can  be 
represented  by  the  opposing  residual  stresses  of  two  continuous  substructures. 

Amenability  of  the  three  approaches  to  multi-axial  phenomena  remains  a  subject  of  investigation.  Promise  in 
overcoming  the  complications  of  mathematically  describing  multi-dimensional  SME  has  been  demonstrated  by 
Liang  and  Rogers,  Boyd  and  Lagoudas,  a  hybrid  micro-plane  approach  proposed  by  Brocca  et  al.,  and  multi¬ 
dimensional  finite  variant  extension  from  the  one-dimensional  Liang  and  Rogers  model  developed  by  Pollard.^"' 

C.  Objectives 

The  objectives  of  this  research  were  to  develop  and  refine  an  SMA  modeling  approach,  the  MTR  subroutine,  to 
predict  the  transient  response  of  a  spacebome  membrane  optic  SMA  deployment  actuator  spine  system  (Fig.  1). 
Motivation  for  this  figure  acquisition  concept  can  be  sought  by  studying  thermonasty,  or  plant  movement  response 
to  temperature.  Tensioning  a  compliant  concentrator  with  thermally  activated  spines  falls  under  the  same  rubric  as 
the  blooming  of  a  flower. 


_  _  ■  concentrator 


Fig.  1  Spaceborne  membrane  optic  SMA  actuator  spine  system  deployment  sequence  concept,  (a)  stowed 
configuration;  (b)  partially  deployed;  (c)  deployed  architecture 

One-time  actuation  of  said  monolithic  actuator  is  anticipated  achievable  by  tailoring  a  SMA  through  variation  of 
elemental  fractions  and  introduction  of  impurities.  The  scope  of  this  study  was  to  model  the  relation  of  input  power 
magnitude  and  waveform  to  stress  fields,  reaction  forces,  and  thermal  fields  for  the  figure  acquisition  of  a  gossamer 
reflector.  Of  particular  interest  is  the  stress  and  thermal  field  history  of  the  polymeric  membrane  concentrator 
through  the  deployment  to  prevent  mechanical  and  thermal  failure  as  limited  experimental  or  modeling  analysis 
results  exist  for  SMA  deployment  schemes. 

This  article  presents  a  unified  approach  from  materials  science,  solid  mechanics,  thermodynamics,  structural 
dynamics,  and  controls  perspectives.  Establishing  accurate  and  practical  methods  for  modeling  the  deployment 
dynamics  of  gossamer  structures,  specifically  to  forecast  spacecraft  loads,  permits  examination  of  different  material 
choices,  and  packaging  strategies,  are  considered  key  to  expedite  the  end  use  of  this  enabling  technology.^ 

II.  MTR  Subroutine 

The  developed  MTR  subroutine  implements  a  phenomenological  constitutive  law  to  model  the  SME  exhibited 
by  SMAs.  The  application  was  initially  designed  to  interface  with  the  ABAQUS  Explicit  solver.  This  finite 
element  tool  supports  a  dynamic,  fully-coupled  thermal-stress  analysis  and  includes  inertia  effects.  The  transient 
thermal  response  is  modeled  by  employing  continuum  elements. Predictions  of  the  shape  memory  effect  require 
the  temperature  and  displacement  degrees  of  freedom  (DOFs)  to  be  solved  simultaneously. 

Acceleration  solutions  are  calculated  at  each  time  step  to  satisfy  dynamic  equilibrium  from  the  equations  of 
motion.  The  kinematic  state  is  advanced  by  employing  a  conditionally  stable,  central  difference  explicit  integration 
rale  and  structuring  the  mass  matrices  diagonally  for  computational  efficiency.  Element  by  element  automatic  time 
incrementation  ensures  the  integration  time  step  is  approximately  equivalent  to  the  least  elemental  dilation  wave 
transient  time  speed.  The  Explicit  package  also  considers  non-linear  geometry,  computes  in  double  machine 
precision,  and  supports  Rayleigh  structural  and  fluid  damping.^^ 

A.  Theory 

Calculation  of  the  martensite  phase  fraction  internal  state  variable,  f,  is  based  on  an  empirical  model  developed 
by  Liang  and  Rogers.^  This  model  accounts  for  the  chemical  free  energy  of  temperature  and  stress,  fitting  the  phase 
transformation  kinetics  to  a  cosine  function  based  on  the  austenite  start  and  finish  temperatures,  and  Af,  and  the 
current  temperature  of  the  alloy,  T;  the  shape  of  the  hysteresis  envelope  resembles  a  sinusoid.  The  form  of  this 
model  implemented  in  the  MTR  subroutine  accounts  for  stress  induced  martensite  (SIM)  through  the  term  driven  by 
the  distortion  energy  stress,  (T„m,  non-dimensionalized  by  the  yield  stress,  Sy,  at  a  rate  inversely  proportional  to  the 
tangent  of  P,  inherently  discounting  superelasticity  occurring  above  the  austenite  finish  temperature  (Eq.  1).  The 
maximum  martensite  phase  fraction  value  of  unity  corresponds  to  complete  rhombohedral  (R)  phase  transformation. 
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The  theory  founding  the  MTR  subroutine  hypothesizes  any  global  geometry  not  exeeeding  the  maximum 
reeoverable  strain  limit  is  a  potential  stable  minimum  energy  state  onee  the  alloy  has  assumed  R  atomie  ordering. 
Conventional  plastie  deformation  is  not  supported.  As  the  SMA  is  subjeeted  to  stress  in  order  to  aehieve  an  altered 
eonfiguration  or  deformed  R  phase  geometry,  the  lattiee  responds  by  reorienting  the  erystalline  strueture  through 
detwinning  and  a  new  global  geometry  stable  minimum  is  sought  eorresponding  to  the  altered  state  of  the  erystal 
lattiee  for  a  given  temperature  and  stress.  For  a  one-way  phase  transformation  from  R  phase  to  body-eentered  eubie 
(B2)  phase,  the  initial  eonfiguration  of  the  finite  element  nodes  modeling  the  SMA  is  eonsistent  with  the  deformed  R 
phase  geometry  stable  minimum.  A  two-phase  erystalline  strueture,  both  R  and  B2,  is  transient  and  global 
geometries  through  the  R  phase  to  B2  phase  transformation  are  eonsidered  unstable. 

It  is  possible  to  relate  the  implementation  of  the  MTR  subroutine  model  to  the  unified  eonstitutive  relation 
published  by  Liang  and  Rogers^  amongst  others,  where  the  stress,  a-og,  depends  on  the  elastie  stress  term,  D{e-Eo,^, 
refereneing  the  eurrent  stable  or  unstable  minimum  geometry,  the  inelastie  stress  term,  Q(f-fo),  aeeounting  for  the 
diserepaney  from  the  trained  shape  or  B2  stable  minimum  geometry  to  the  current  stable  or  unstable  minimum 
geometry,  and  the  thermoelastic  stress  term,  0{T-Tg,^,  referencing  the  current  stable  or  unstable  minimum  geometry 
(Eq.  2). 


<t-<t,=D(e-e„^)  +  C1(P  -  4)  +  0(r  -  To,  (2) 

The  tensors  can  represent  one-,  two-,  or  three-dimensional  stress  fields.  If  the  continuum  is  initialized  at  a 
mechanically  and  thermally  unloaded  equilibrium,  R  phase  configuration  for  a  B2  phase  stable  minimum  geometry 
recovery  analysis,  both  second  order  tensor  naught  values  would  be  zero,  (fo=l,  and  To  would  be  less  than  A^.  The 
elastic  modulus  as  a  function  of  the  martensite  phase  fraction  internal  state  variable  is  defined  as  D{e-Eo,^.  This 
term  is  modeled  continuously  through  the  ABAQUS  linear  elastic  material  law. 

The  detwinned  stress,  also  a  function  of  the  martensite  phase  fraction  internal  state  variable,  is  defined  as 
and  is  modeled  discreetly  by  updating  the  nodal  geometry  vectors  between  successive  loop  analyses.  As  the  model 
retransforms  back  to  the  B2  phase,  the  nodal  geometry  vectors  mapping  neighboring  nodes  track  from  the 
configuration  consistent  with  the  R  phase  global  stable  minimum. 


where 


i  =  l,2,...,A„-l, 

to  the  configuration  consistent  with  the  B2  phase  global  stable  minimum. 


where 


i  =  \X-,N^-\, 


forming  a  set  of  unstable  minimum  energy  states  through  the  transformation  along  paths  eoineident  with  both  R  and 
B2  phase  geometries.  Mimieking  MTR  by  defining  these  transformation  veetors  is  requisite  for  a  high  fidelity 
model  as  there  is  only  one  atomistie  path  from  lower  to  higher  symmetry  between  the  R  and  B2  phase  erystalline 
struetures.  This  term,  Q(f-fo),  in  effeet,  realizes  and  imposes  inereased  detwinned  strain  as  the  martensite  phase 
fraetion  internal  state  variable  deereases, 

-G,„)+(1-^,,,,Xg,,,,  -G,J,  (3) 

where 

and 

j  =  l2,...,N,, 

and  eonsequently  drives  the  return  to  the  parent  eonfiguration  with  no  assoeiated  detwinned  stress,  i.e.,  the  trained 
austenite  shape. 

This  eonstituent  hypothesis  is  the  erux  of  the  subroutine’s  approaeh  to  emulate  the  SME.  Mapping  originates 
from  nodes  of  known  eonfiguration;  reviewing  Eq.  3,  node  i-l  would  be  deelared  at  an  essential  boundary. 
Progress  along  the  transformation  paths  is  proportional  to  the  martensite  phase  fraetion  internal  state  variable.  This 
set  of  unstable  minima  is  not  strietly  globally  deterministie  sinee  the  paths  the  nodal  geometry  veetors  traek  are 
relative  to  the  eonfiguration  of  neighboring  nodes,  but  these  veetors  do  have  a  strong  attraetion  to  the  B2  phase 
global  stable  minimum  geometry.  Elemental  diseretization  of  the  martensite  phase  fraetion  supports  unique 
transformation  paths  and  eonsequently  diseretization  of  the  SME  to  the  finite  element  seale.  The  relative  nature  of 
the  nodal  geometry  veetors  does  not  eompromise  global  deployment  destiny  sinee  only  one  parent  phase  global 
stable  minimum  eonfiguration  is  possible. 

The  MTR  subroutine  eapability  is  limited  to  one-way  phase  transformations  from  martensite  to  austenite;  two- 
way  and  superelastieity  behavior  is  not  supported.  As  the  subroutine  was  designed  to  be  implemented  in 
eonjunetion  with  a  deployment  model  not  to  exeeed  the  threshold  stress,  these  phenomena  were  not  of  eoneem. 
Consistent  with  the  proposed  theory,  the  subroutine  reassigns  nodal  geometry  veetors,  as  well  as  phase  dependent 
properties,  ineluding  modulus  of  elastieity,  CTE,  and  thermal  eonduetivity,  proportional  to  the  Liang  and  Rogers 
martensite  phase  fraetion  internal  state  variable  between  ABAQUS  analyses  (Eqs.  4-6).  The  duration  of  these 
analyses  is  orders  of  magnitude  less  than  the  period  of  the  model’s  first  austenitie  mode  to  avoid  infiuenee 
sensitivity  at  the  natural  frequeney  and  ensure  suffieient  temporal  diseretization. 

v  )=  4/ K  (4) 

where 

!=1,2,...,A^ 

j  =  ia,-,N, 

fe-.y  )  =  f,-.y«M  +  (l  -  ^/.y 

where 


i  =  l2,...,N, 


j  =  l2,...,N, 


where 


(6) 


i  =  \,2,-,N^ 
j  =  l2,...,N, 

The  subroutine  maintains  displacements  first  order  continuously  differentiable  and  temperatures  zeroth  order 
continuously  differentiable  by  imposing  boundary  conditions  of  displacement  and  velocity,  and  initial  conditions  of 
temperature,  at  the  beginning  of  each  analysis  transcribed  from  the  final  values  of  the  previous  analysis.  Boundary 
condition  steps  are  required  for  the  mechanical  DOFs  because  the  ABAQUS  Explicit  solver  will  not  permit 
displacement  discontinuities.  These  conditions  are  imposed  over  durations  many  orders  of  magnitude  less  than  the 
period  of  the  model’s  first  austenitic  mode  to  minimize  their  influence  on  the  simulation. 

B.  Architecture 

As  discussed,  spatial  discretization  of  the  SME  is  to  the  elemental  level  of  the  finite  element  model  (FEM)  and 
temporal  discretization  is  to  the  interval  of  the  ABAQUS  analyses.  Recall  the  subroutine  updates  the  relative  nodal 
geometry  vectors  discretely  between  ABAQUS  analyses  to  realize  and  impose  increased  detwinned  strain 
proportional  to  the  martensite  phase  fraction  internal  state  variable.  A  simulation  run  is  composed  of  many 
ABAQUS  analyses  and  MTR  subroutine  executions  as  illustrated  in  the  loop  map  (Fig.  2). 
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Fig.  2  ABAQUS  aud  MTR  subroutiue  simulatiou  ruu  loop 
architecture 


A  simulation  run  for  a  given  model  is  primed  with  an  ABAQUS  analysis  of  an  input  file  enforcing  the  initial 
conditions.  Typically,  the  SMA  components  of  the  FEM  are  assumed  to  be  completely  martensitic  at  the 
commencement  of  the  run;  the  initial  configuration  of  their  respective  finite  element  nodes  is  consistent  with  the 
deformed  R  phase  geometry  stable  minimum.  A  unique  element  set  and  material  property  definition  is  associated 
with  each  SMA  element  to  allow  discretization  of  the  SME  to  the  finite  element  scale  for  the  four  material  properties 
dependent  upon  the  phase  transformation.  ABAQUS  Explicit  permits  temperature  initial  conditions  to  be  enforced. 
These  DOFs,  generally  initialized  below  the  B2  phase  start  temperature,  are  imposed  on  all  the  nodes  of  the  model 


before  the  solver  analyzes  the  steps.  This  component  of  the  input  file  maintains  temperatures  zeroth  order 
continuously  differentiable  between  subsequent  loop  analyses. 

The  first  step  is  composed  of  a  displacement  boundary  condition  with  authority  over  all  the  nodes  of  the  model. 
This  component  of  the  input  fde  enforces  any  initial  elastic  strain  and  maintains  displacements  zeroth  order 
continuously  differentiable  between  subsequent  loop  analyses.  The  second  step  of  the  input  file  is  composed  of  a 
velocity  boundary  condition  with  authority  over  all  the  nodes  of  the  model.  This  component  of  the  input  file 
enforces  any  initial  strain  rate  and  maintains  displacements  first  order  continuously  differentiable  between 
subsequent  loop  analyses. 

The  magnitude  of  both  these  essential  boundary  conditions  is  proportional  to  a  linear  amplitude  profile  defined 
by  a  function  coincident  with  zero  and  unity  at  the  beginning  and  end  of  the  step,  respectively.  To  optimize  for 
computational  efficiency,  the  duration  of  the  first  two  steps  is  generally  set  at  one  order  of  magnitude  greater  than 
the  least  elemental  dilation  transient  time  wave  speed  of  the  FEM.  Modeling  activity  at  rates  greater  than  this  FEM 
property  is  not  communicated  through  the  continuum  to  be  included  in  the  equilibrium  calculations  and  results  in  an 
unstable  solution.  Both  of  these  steps  respect  the  mechanical  essential  boundary  conditions  defined  for  the  model. 

The  third  step  is  the  only  step  of  an  individual  ABAQUS  analysis  contributing  to  the  simulation  run  solution.  In 
effect,  each  third  step  solution  can  be  considered  a  finite  segment  of  the  model  history.  In  addition  to  also 
respecting  the  mechanical  essential  boundary  conditions,  thermal  essential  and  mechanical  and  thermal  natural 
boundary  conditions,  including  tractions  and  body  loads  as  defined  for  the  model,  are  included.  The  duration  of  this 
step  is  often  found  between  competing  trends  of  SME  temporal  discretization  demands  and  computational  capacity. 

Upon  completion  of  an  ABAQUS  analysis,  the  binary  results  file  is  converted  to  ASCII  format  to  be  read  by  the 
MTR  subroutine.  The  first  component  of  the  subroutine  captures  the  nodal  configurations.  Through  practice  it  was 
discovered  ABAQUS  does  not  transcribe  nodal  geometry  vectors  from  the  input  file  to  the  results  file  in  double 
precision.  To  maintain  simulation  run  double  precision,  the  nodal  declarations  are  read  directly  from  the  input  file. 

The  second  component  of  the  MTR  subroutine  captures  final  third  step  nodal  displacements,  velocities, 
accelerations,  temperatures,  and  reaction  forces  and  elemental  distortion  energy  stresses  from  the  converted 
ABAQUS  results  file.  After  a  set  interval  of  simulation  run  loops,  Ni,  dictated  by  model  history  resolution 
requirements,  these  solutions  are  appended  to  the  simulation  run  solution  log  results  file.  The  logged  displacements 
are  manipulated  to  yield  identical  nodal  positions  as  displaced  from  the  B2  phase  global  stable  minimum  geometry; 
this  vector  is  defined  as  the  mapping  of  the  configuration  resulting  from  inelastic  strain  and  displacements  from 
elastic  strain. 


from  the  B2  phase  geometry. 


(7) 


where 


1=1,2,.. 
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n 


and 


j  =  l2,...,Nr 

Both  R  and  B2  phase  nodal  configurations  are  embedded  in  the  subroutine.  A  log  file,  with  displacement  results 
applicable  to  a  single  model,  enables  continuous  post  processing  of  the  simulation  run  solution  as  opposed  to 
processing  the  ABAQUS  analyses  third  step  results  individually,  with  the  respective  unique  nodal  declarations. 

To  prepare  for  calculating  the  Liang  and  Rogers  martensite  phase  fraction  internal  state  variable  for  each 
element,  the  subroutine  distills  an  average  elemental  temperature  from  the  temperatures  of  the  nodes  defining  the 
element.  B2  phase  start  and  finish  temperatures,  B2  phase  yield  stress,  and  the  SIM  angle,  as  well  as  elemental 
initial  martensite  phase  fraction  are  embedded  in  the  subroutine.  With  all  the  R  phase  fraction  variables  accounted, 
Eq.  1  is  solved  for  each  element  followed  by  Eq.  3  for  each  node.  As  discussed,  the  relative  nodal  geometry  vector 
formulation  respects  the  mechanical  constraints  defined  for  the  model.  With  this  approach,  geometries  of 


constrained  nodes  are  considered  constants,  serving  as  a  ground  for  the  calculation  of  the  unconstrained  node 
declarations.  Phase  dependent  properties  are  then  solved  for  at  the  finite  elemental  scale  based  on  the  R  and  B2 
phase  values  also  embedded  in  the  subroutine  as  Eqs.  4-7. 

Upon  executing  these  tasks,  the  next  component  of  the  MTR  subroutine  writes  an  ABAQUS  input  fde  with  the 
three  step  architecture  to  enter  into  the  next  loop  of  the  simulation  run.  Loops  are  completed  until  the  model  history 
of  interest  is  solved. 

C.  Validation 

With  the  MTR  subroutine  operational,  experimental  studies  of  beam  test  articles  were  conducted  to  validate  the 
fidelity  of  the  MTR  subroutine  and  support  tuning  of  parameters  governing  structural  damping  and  the  SIM 
phenomenon.  These  deflected  R  phase  beam  specimens  were  subjected  to  a  series  of  temperature  ramp  magnitudes 
to  represent  various  body  flux  or  ohmic  heating  thermal  inputs.  Restraining  force  loads  were  applied  at  the  deflected 
boundary  to  invoke  responses  under  various  expected  in-service  functional  specifications. 

1.  Experimental  Design 

A  TA  Instruments  Q800  dynamic  mechanical  analyzer  (DMA)  was  configured  with  a  single  cantilever  clamp 
fixture  for  the  experimental  studies  to  simulate  a  scaled  membrane  optic  spine  deployment  actuator.  The  DMA 
single  cantilever  clamp  arrangement  was  chosen  as  the  most  representative  of  the  envisioned  figure  acquisition 
reproduced  in  Fig.  1.  This  instrument  was  specified  at  a  force  resolution  of  100  pN,  a  displacement  resolution  of 
1 .00  nm,  and  an  isothermal  stability  of  100  mK.^^ 

The  model  considered  to  validate  the  fidelity  of  the  MTR  subroutine  and  support  tuning  of  parameters  governing 
structural  damping  and  the  SIM  phenomenon  was  a  fixed-roller  boundary  beam.  This  two-dimensional  system 
allows  for  realizing  relatively  high  strain  potential  to  generated  force  and  is  applicable  to  many  SMA  actuation 
schemes.  In  addition,  this  setup  is  favorable  to  address  life-cycle  design  constraints  of  SMA  actuators,  namely  the 
strain  to  fatigue  degradation  relation.  The  dimensional  scale  of  the  test  specimen  is  a  gauge  length  of  17.4  mm  (Fig. 
3). 


y 


Fig.  3  Fixed-roller  boundary  beam  actuator 
model 


The  SMA  selected  for  experimental  work  on  the  basis  of  square  cross  section  availability  was  SE-508,  a 
commercially  available  Nitinol  alloy.  A  square  cross  section  allowed  for  an  initial  coarse  mesh  representation  of  the 
specimen  with  brick  elements.  SE-508  strip  with  a  660  pm  square  cross  section  and  black  oxide  surface  finish  were 
procured. 

The  specimen  was  subjected  to  a  series  of  16.7  mKs  ’  and  167  mKs  '  linear  temperature  ramp  magnitudes  to 
simulate  varying  body  flux  or  ohmic  heating  thermal  loads,  Qbj(r,t).  Magnitudes  greater  than  167  mKs  *  resulted  in 
pronounced  non-uniform  heating  of  the  specimen.  In  addition  to  the  free  recovery  scenario,  constant  1.00  N  and 
10.0  N  restraining  force  loads,  F,  were  applied  at  the  roller  boundary  parallel  to  the  y-axis  impeding  the  global 
geometry  recovery  to  simulate  return  with  expected  in-service  functional  specifications,  such  as  the  inertial  load  of 
the  reflector,  as  well  as  other  imaging  hardware. 


The  procedure  followed  for  each  scenario  began  by  mounting  the  beam  test  articles  in  the  clamps  torqued  to  1.10 
Nm.  The  moveable  clamp  was  then  freed  and  floated  to  allow  the  specimen  to  assume  its  self-weight  shape.  The 
moveable  clamp,  previously  zeroed  during  the  calibration  sequence,  was  locked.  This  step  permitted  the 
examination  of  gravity  accelerations  on  the  periphery  displacement.  Effects  due  to  this  body  load  were  found  to  be 
approximately  116  nm  and  considered  negligible. 

The  furnace  was  then  cooled  below  the  R  phase  finish  temperature  and  maintained  at  this  themial  condition  for 
180  s.  After  the  delay,  a  10.0  N  load  was  engaged  at  the  roller  boundary  in  the  direction  of  the  gravity  vector  or  the 
negative  y-unit  vector  and  paused  for  180  s.  This  load  was  then  either  disengaged  or  adjusted  to  the  magnitude  of 
the  restraining  force  and  maintained  for  180  s.  This  periphery  displacement  was  assumed  to  be  characteristic  of  the 
deformed  R  phase  geometry  of  the  beam.  (Note:  the  detwinned  configuration  for  all  scenarios  qualified  within  the 
range  of  recoverable  inelastic  strains  for  Nitnol.®) 

After  this  third  delay  the  linear  temperature  ramp  load  was  engaged  and,  if  applicable,  the  restraining  force  was 
maintained.  Recorded  temperature  and  displacement  history  from  this  segment  of  the  test  procedure  was  to  be 
modeled  in  the  simulation  run  to  support  tuning  of  parameters  governing  structural  damping  and  the  SIM 
phenomenon.  Plots  of  roller  boundary  displacement  history  and  roller  boundary  force  against  displacement  for  a 
scenario  characterized  by  a  167  mKs  *  temperature  ramp  and  1.00  N  restraining  force  are  illustrative  of  the  test 
procedure  (Figs.  4,  5). 
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Fig.  4  Roller  boundary  displacement  history  experimental  results  of  167  mKs  *  temperature  ramp  and  1.00  N 

restraining  force  scenario 
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Fig.  5  Roller  boundary  force  against  displacement  experimental  results  of  167  mKs  *  temperature  ramp  and 

1.00  N  restraining  force  scenario 


2.  Model  Refinement 

The  MTR  subroutine  was  customized  to  the  fixed-roller  boundary  beam  model  to  replicate  the  DMA 
experiments  by  varying  expected  influential,  unknown  parameters  governing  the  SIM  phenomenon  and  structural 
damping.  The  FEM  was  composed  of  108  nodes.  A  three-dimensional  model  was  considered  to  support 
development  of  the  subroutine’s  utility  for  multi-dimensional  behavior.  The  envisioned  figure  acquisition 
reproduced  in  Fig.  1  would  involve  radial  spines  and  not  a  single  axisymmetric  actuator.  The  108  nodes  defined  26 
ABAQUS  C3D8RT  elements  as  close  to  equal  sided  cubes  as  single  element  depth  model  representation  would 
allow.  Recognizably,  the  model  was  coarsely  meshed  to  minimize  computational  expense.  The  B2  phase  geometry 
is  defined  by  Fig.  3.  Material  properties  are  consistent  with  those  determined  experimentally  or  published  values. 

27 

At  the  onset  it  was  apparent  it  would  not  be  feasible  to  employ  the  ABAQUS  Explicit  solver  for  simulations  of 
this  duration  due  to  the  computational  demanding  integration  method.  The  MTR  subroutine  was  modified  to 
interface  with  the  ABAQUS  Standard  solver.  Direct  integration  solutions  are  calculated  at  each  time  step  to  satisfy 
dynamic  equilibrium.  The  kinematic  state  is  advanced  by  employing  the  implicit  Hilber-Hughes-Taylor  operator 
developed  from  the  trapezoidal  rule,  implemented  iteratively  with  Newton’s  method  for  convergence  of  a  non-linear 
equation  set.^*  Implications  of  this  retooling  to  the  thermal  solution  and  timing  aspects  of  the  simulation  run  are 
discussed  in  Ref  24. 

These  simulation  runs  were  solved  on  a  32.0  bit,  5.60  GHz  dual  processor  machine  with  2.00  GB  of  RAM. 
Approximately  120  machine  hours  were  required  to  complete  a  run.  This  investment  could  have  been  reduced 
approximately  60.0  %  by  representing  the  test  article  with  a  two-dimensional  model  if  three-dimensional 
performance  was  not  of  concern. 

The  scenario  characterized  by  no  restraining  force  was  the  first  attempted  to  replicate.  SIM  behavior  was 
discounted  by  valuing  P  at  1.57  radians.  Rayleigh  stiffness  proportional  damping  was  enforced  at  a  nominal  value 
of  10.0  %  of  critical  relative  to  the  first  austenitic  mode. 

In  order  to  improve  model  fidelity,  the  Liang  and  Rogers  phase  transformation  kinetics  model  was  modified  with 
an  additional  series  of  seven  correction  terms  to  capture  the  features  of  the  SME  distribution  exhibited  by  the  test 
article.  Although  the  single  term  formulation  of  Eq.  1  may  mimic  the  SME  distribution  of  a  single  crystal  SMA,  this 
idealization  can  not  be  expected  of  a  polycrystalline  SMA  actuator  in  practice.  These  sinusoidal  correction  terms, 
with  parameters  found  from  a  least  squares  fit  to  the  experimental  results  to  seek  a  more  accurate  simulation  run 
solution,  were  typically  cast  in  pairs  to  ensure  continuity  at  the  engagement  and  disengagement  of  the  Liang  and 
Rogers  term.^"'  The  martensite  phase  fraction  correction  terms  were  incorporated  into  the  MTR  subroutine  and  the 
replication  was  again  executed  (Fig.  6). 


temperature, 

K 


Fig.  6  MTR  subroutine  solution  of  roller  boundary  displacement  bistory  of  167  mKs  *  temperature  ramp  and 
no  restraining  force  scenario  with  martensite  phase  fraction  internal  state  variable  correction  terms 


As  the  stress  state  of  the  specimen  was  not  significant,  little  could  be  gained  from  this  scenario  regarding  the 
SIM  parameter.  With  regard  to  structural  damping,  dependencies  on  this  energy  dissipation  term  were  investigated 
more  thoroughly  considering  the  restraining  force  scenarios.  The  value  in  this  scenario  was  the  reformulation  of  the 
Liang  and  Rogers  phase  transformation  kinetics  model.  With  this  refinement  coded,  efforts  focused  on  replicating 
the  restraining  force  experiments.  (Note:  all  subsequent  simulation  runs  incorporated  the  seven  martensite  phase 
fraction  correction  terms  in  addition  to  the  original  Liang  and  Rogers  phase  transformation  kinetics  model  into  the 
MTR  subroutine.) 

The  scenario  characterized  by  a  1.00  N  restraining  force  was  the  next  attempted  to  replicate.  SIM  behavior  was 
initially  discounted  by  valuing  P  at  1.57  radians.  Again,  Rayleigh  stiffness  proportional  damping  was  enforced  at  a 
nominal  value  of  10.0  %  of  critical  relative  to  the  first  austenitic  mode.  This  simulation  run  encountered  a  numerical 
instability  at  approximately  100  s  simulation  time  and  faulted  as  the  elements  declarations  were  too  distorted  for 
convergence  before  any  meaningful  insight  could  be  obtained. 

Reference  24  documents  an  investigation  of  the  numerical  artifacts  which  arise  when  modeling  constant  natural 
boundary  conditions  in  attempt  to  map  the  SIM  parameter  window.  Research  performed  by  Bo  et  al.^^  proved 
insightful  as  to  the  experimentally  measured  value  of  this  Nitinol  material  property.  Extrapolating  from  the 
presented  stress  against  temperature  plots,  the  value  of  the  SIM  angle  was  interpreted  as  approximately  1 .47  radians. 
In  order  to  compare  with  the  values  coded  in  the  MTR  subroutine,  this  angle  must  be  refigured  for  a  non- 
dimensional  relationship.  The  non-dimensional  empirical  value  of  the  angle  is  approximately  1.41  radians  with  358 
K  and  560  MPa  as  the  values  of  the  B2  phase  finish  temperature  and  the  B2  phase  yield  stress,  respectively. 
Reference  24  also  assesses  the  influence  of  structural  damping;  the  MTR  subroutine  solutions  were  found  to  be 
insensitive  to  this  parameter  at  the  experimental  phase  transformation  rates. 

A  zero  gravity  variant  of  the  scenario  characterized  by  a  1 .00  N  restraining  force  was  executed.  The  constant 
boundary  condition  at  the  periphery  was  replaced  with  an  equivalent  mass  of  the  1.00  N  weight  force  within  a 
terrestrial  gravitational  field  or  0.102  Kg.  It  was  anticipated  this  inertial  loading  scenario  would  be  representative  of 
the  deployment  sequence  predictions  and  consequently,  it  was  necessary  to  determine  any  attributable  instabilities  or 
divergence  trends.  SIM  behavior  was  accounted  for  by  valuing  P  at  1.37  radians.  Again,  Rayleigh  stiffness 
proportional  damping  was  enforced  at  a  nominal  value  of  10.0  %  of  critical  relative  to  the  first  austenitic  mode  (Fig. 
7).  (Note:  this  simulation  run  yielded  an  intuitive  roller  boundary  displacement  history,  with  any  numerical 
complications  noticeably  absent.) 
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Fig.  7  MTR  subroutine  solution  of  roller  boundary  displacement  bistory  of  167  mKs  *  temperature  ramp  and 

0.102  Kg  mass  scenario 


III.  Deployment  Sequence 


A.  Deployment  Model  Design 

The  MTR  subroutine  was  eustomized  to  represent  an  individual  SMA  radial  aetuator  and  eorresponding 
eoneentrator  sector  of  a  spaceborne  membrane  optic  SMA  deployment  actuator  spine  system  to  forecast  the 
envisioned  figure  acquisition  reproduced  in  Fig.  1  (Fig.  8).  The  imagined  system  architecture  is  composed  of  36 
circumferentially  equally  spaced  SMA  radial  actuators.  This  model  has  the  parent  geometry  of  a  parabolic  profiled 
spine  parallel  to  the  x-axis  with  the  nominal  dimensions  of  1.00  m  projected  area  radius,  Rs,  10.0  mm  square  cross 
section  thickness,  hs,  and  a  focal  length  to  diameter  ratio  /x.,s=0.500  peripherally  supporting  a  parabolic,  polymeric 
reflector  with  the  nominal  dimensions  of  1.00  m  projected  area  radius,  Rc,  25.0  pm  thickness,  he,  175  mm  width, 
Wq,  and  a  focal  length  to  diameter  ratio  equal  to  unity.  This  shape  is  considered  to  be  the  equilibrium 
configuration  of  the  aperture  continuum  reflecting  current  manufacturing  concepts.^** 


Rc—Rs 


Fig.  8  Deployment  spine  and  concentrator 
model 


The  cross  section  of  the  concentrator  sector  was  considered  prismatic  along  the  radius  or  parallel  to  the  x-axis  for 
this  study  and  the  curvature  of  the  periphery  boundary  was  approximated  by  the  tangent  to  the  polymeric  reflector 
projected  area  at  the  spine  interface.  Rigid,  beam,  three-dimensional,  and  two-node  (ABAQUS  RB3D2)  elements 
were  modeled  at  this  boundary  appended  to  the  SMA  spine,  effectively  serving  as  the  mechanical  coupling  between 
the  two  continuums;  circumferential  load  paths  were  ignored  by  modeling  the  rigid  elements  along  the  tangent. 
These  model  reductions  do  not  prohibit  investigation  of  the  transient  thermal  and  radial  stress  states  of  the 
concentrator,  expected  to  be  of  greatest  concern. 

The  FEMs  were  composed  of  the  legacy  108  nodes  defining  26  ABAQUS  C3D8RT  elements  format 
representing  the  spine  and  162  nodes  defining  130  shell,  four-node,  and  reduced  integration  (ABAQUS  S4R) 
elements  representing  the  concentrator  (Fig.  9).  A  three-dimensional  model  was  considered  as  the  deployment 
involves  multi-dimensional  behavior.  Envisioned  figure  acquisition  reproduced  in  Fig.  1  would  involve  radial 
spines  and  not  a  single  axisymmetric  actuator.  Recognizably,  the  model  was  coarsely  meshed  to  minimize 
computational  expenses. 


Fig.  9  FEM  representation  of  the  deployment  spine  and  concentrator  model  with  rigid  elements  at 
concentrator  periphery  honndary  appended  to  the  SMA  spine 

Spine  material  properties,  ^5=33. 1-57.0  GPa,  U5=0.330,  as=6.62xl0'^-1.20xl0'^  K"',  p5=6,500  kgm'^,  and 
^^5=8. 60-18.0  Wm'*K'*,  are  of  SE-508,  a  representative  Nitinol  alloy.  Concentrator  material  properties,  iic^2.80 
GPa,  U(-=0.340,  ac^2.00xl0"^  K  ',  pc^l,420  kgm'^  Cc=l,090  Jkg  'K  ',  and  i:c^0.120  Wm  'K  ',  are  of  Kapton  100 
HN®,  representative  of  a  polyimide  material.  Phase  dependent  properties  independently  sought  include  spine 
modulus  of  elasticity,  E^,  and  CTE,  a*  by  employing  a  DMA  and  TMA  instruments,  respectively.  Thermal 
conductivity  was  assumed  to  be  bounded  with  published  values. In  order  to  accurately  predict  the  thermal- 
mechanical  response  of  any  specimen  of  this  alloy,  it  is  requisite  to  determine  these  properties,  as  well  as  the  B2 
phase  start  and  finish  temperatures,  from  DSC  testing  of  the  specific  material  in  hand.^"^ 

In  order  to  model  thermal  loads,  effective  specific  heat  as  a  function  of  temperature  must  be  invoked  for 
continuums  undergoing  phase  changes.  The  term  effective  specific  heat  implies  both  the  latent  energy  of  the 
martensite  to  austenite  phase  transformation  and  the  sensible  heat  rate  are  accounted  for. 

In  addition  to  B2  start  and  finish  temperatures,  effective  specific  heat,  Cs,  was  distilled  from  SE-508  DSC 
results.  This  material  property  as  a  function  of  temperature  was  found  by  numerically  integrating  the  specific  heat 
flow  to  gain  a  function  of  specific  energy  addition  as  a  function  of  temperature.  Numerically  differentiating  this 
intermediate  solution  by  temperature  results  in  the  change  in  specific  energy  addition  against  temperature  or 
effective  specific  heat.  This  relationship  was  fit  to  a  second  degree  polynomial  to  be  embedded  in  the  MTR 
subroutine  (Eq.  8). 


c,(r)  =  -o.o288r"^!— +  i8.9r^!— -2340^!— 

KgK^  KgK'  KgK 

Natural  and  essential  boundary  conditions,  Qbjir,t),  Qc{r,t),  P,  and  plane  symmetry  of  each  model  continuum, 
were  enforced  to  simulate  realistic  deployment  scenarios.  Three  specific  thermal  body  flux  load  waveforms,  qi,/r,t), 
were  considered  within  current  projected  on-orbit  power  budgets  corresponding  to  three  loading  scenarios  (Eqs.  9- 
11).^'  Thermal  constraints  restricted  convective  fluxes  at  all  surfaces,  i.e.,  Qc(r,t)^0,  and  radiation  effects  were 
ignored. 
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Consistent  with  general  eontrol  methodology  of  membrane  opties,  all  seenarios  ineluded  a  uniform,  bias,  1.00 
Pa,  surfaee-normal,  and  eonstant  load  aeting  on  the  coneentrator.  This  load  eould  be  provided  by  gas  or,  to  avoid 
attendant  makeup  gas  and  lentieular  eomplieations,  eleetrostatie  pressure.  It  is  expeeted  some  degree  of  self- 
rigidization  will  be  required  to  not  only  aequire  the  deterministie  figure,  but  also  for  maintenanee  of  the  polymerie 
refleetor’s  figure,  even  with  the  aid  of  a  strain  eoating.^**  The  aperture  was  meehanieally  eonstrained  from 
translations  in  the  direetions  of  the  x-  and  z-unit  veetors,  with  translations  in  the  direetion  of  the  y-unit  veetor  and 
rotational  DOFs  uneonstrained,  at  the  vertex  boundary  and  the  spine  was  meehanieally  eonstrained  from  all  three 
translations  at  the  vertex  boundary,  eonsistent  with  the  quasi-axisymmetrie  reduetion  of  the  model. 

The  R  phase  geometry  of  the  SMA  spine  and  deformed  state  of  the  polymerie  refleetor  were  found  by  elastieally 
translating  the  two  interfaee  nodes  at  the  periphery,  (1.000  m,  0.260  m),  of  the  model  in  the  direetion  of  the  negative 
x-unit  veetor  to  the  vertieal  axis,  (0.000  m,  0.81 1  m),  while  not  permitting  translations  of  these  nodes  parallel  to  the 
z-axis  (Fig.  10).  Shell,  four-node,  and  full  integration  (ABAQUS  S4)  elements  were  utilized  and  eontaet 
phenomenon  was  diseounted.  The  ABAQUS  statie  solver  was  employed.  This  finite  element  analysis  maintained 
the  essential  boundaries  at  the  respeetive  vertiees  of  eaeh  eontinuum  and  the  natural  boundary  eondition  preseribed 
at  the  surfaee  of  the  eoneentrator,  with  non-linear  deformation  invoked.  This  global  geometry  is  eharaeterized  by  a 
maximum  prineipal  spine  strain  of  8.43  %  oeeurring  in  element  26  at  the  periphery  and  qualifying  within  the 
maximum  reeoverable  strain  limit  for  one  time  aetuation  of  Nitinol.^  Considered  as  the  stowed  eonfiguration,  this 
paekaging  reduees  the  projeeted  area  radius  of  the  observatory  by  approximately  one-third. 


Fig.  10  R  phase  geometry  of  the  SMA  spine  and  deformed  state  of  the  polymeric 

reflector 


Based  on  the  eomputational  eapaeity  limitations  eneountered  with  refinement  efforts,  whieh  precluded  the 
feasibility  of  employing  the  ABAQUS  Explicit  solver  for  simulations  of  the  duration  of  the  experimental  scenarios 
due  to  the  computational  demanding  integration  method,  it  was  anticipated  the  ABAQUS  Standard  solver  would  be 
utilized  for  the  deployment  sequence  predictions  as  well.  Implications  of  this  retooling  to  the  thermal  solution  and 
timing  aspects  of  the  simulation  run  are  discussed  and  an  additional  complexity  associated  with  this  simulation  run 
as  the  model  comprises  of  components  which  do  not  undergo  a  change  in  static  equilibrium  state,  i.e.,  the 
concentrator,  is  addressed  in  Ref.  24. 


SIM  behavior  was  accounted  for  by  valuing  P  at  1.37  radians.  Rayleigh  stiffness  proportional  or  viscous 
material  damping  was  enforced  at  10.0  %  of  critical  relative  to  the  first  modes  of  both  the  austenitic  spine  and 
concentrator.  Enforcing  energy  dissipation  through  this  linear  construct  accounts  for  internal  friction  of  the  original 
crystallographic  structure  recovery  of  the  spine  and  possible  radiative  damping,  a  manifest  of  the  bias  pressure  load 
on  the  concentrator.  Reference  24  details  further  development  of  the  MTR  subroutine  relevant  to  forecasting  the 
deployment  sequence. 

B.  Deployment  Model  Results 

The  scope  of  this  study  was  to  model  the  relation  of  input  power  magnitude  and  waveform  to  stress  fields, 
reaction  forces,  and  thermal  fields  for  the  figure  acquisition  of  a  gossamer,  polymeric  reflector.  Of  particular 
interest  is  the  stress  and  thermal  field  history  of  the  polymeric  membrane  concentrator  through  the  deployment  to 
prevent  mechanical  and  thermal  failure  as  limited  experimental  or  modeling  analysis  results  exist  for  SMA 
deployment  schemes. 

The  MTR  subroutine  simulation  run  solution  enforcing  the  step  waveform  body  flux  loading  scenario  defined  by 
Eq.  9  yielded  an  intuitive  global  response.  Stills  of  the  SMA  spine’s  detwinned  strain  recovery  are  illustrated  at  57.6 
s  and  115  s  simulation  run  time  (Fig.  11).  It  is  evident  from  these  unstable  minimum  energy  state  configurations, 
logic  implementing  compatibility  theory  to  ensure  conservation  of  mass  is  required  for  greater  fidelity  of  SME 
predictions  undergoing  an  R  to  B2  phase  global  geometry  return,  with  relatively  large  attendant  translations  and 
rotations. 


(a) 


Fig.  11  Spine  FEM  detwinning  strain  recovery  responding  to  Eq.  9  step  waveform  body  flnx  loading 
scenario,  (a)  at  57.6  s  simnlation  rnn  time;  (b)  at  115  s  simnlation  rnn  time 


The  figure  acquisition  simulation  run  solutions  enforcing  step  waveform  body  flux  loading  scenarios  periphery 
displacement  histories  present  a  more  detailed  perspective  of  the  deployment  (Fig.  12).  Temperature  history  of  the 
interface  is  also  plotted. 


Fig.  12  MTR  subroutine  solution  of  periphery  x-  and  y-displacement  and  temperature  history  responding 

to  Eq.  9  step  waveform  body  flux  loading  scenario 


The  vibration  behavior  present  in  all  three  y-displacement  histories  preceding  the  preliminary  y-component 
equilibrium  excursions  range  from  0.077  Hz  to  0.160  Hz  for  a  four  cycle  average  of  the  responses  excited  by  the 
load  waveform  scenarios  defined  by  Eqs.  10,  11,  respectively.  The  frequency  of  the  first  austenitic  mode  of  the 
spine  in  the  unloaded  equilibrium  configuration,  divorced  from  the  concentrator,  was  found  to  occur  at  0.521  Hz. 
One  would  expect  the  natural  frequency  of  the  spine  and  concentrator  integrated  model  in  an  intermediate  phase 
state  to  be  less  than  this  property,  with  the  mass  of  the  concentrator  effecting  at  the  periphery  and  devalued  elastic 
modulus.  The  period  of  the  first  excursion  cycle  following  aperture  tensioning  for  all  three  scenarios  was 
approximately  90.0-95.0  s,  most  probably  the  natural  period  of  the  integrated  model  excited  by  the  concentrator 
engagement. 

Distortion  energy  stress  history  is  reported  at  the  periphery  of  the  spine  from  the  three  simulation  run  solutions 
(Fig.  13).  These  contours  are  the  running  sum  of  the  elastic  and  inelastic  stress  terms  of  the  unified  constitutive 
relation  introduced  as  Eq.  2.  (Note:  according  to  published  values  these  stress  spikes  upon  aperture  tensioning  are 
within  Nitinol’s  austenite  proportional  limit  range  to  a  maximum  of  6.90x10*  Pa.^’) 
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Fig.  13  MTR  subroutine  solution  of  spine  periphery  distortion  energy  stress  bistory  responding  to  Eqs.  9- 

11  waveform  body  flux  loading  scenarios 


Distortion  energy  stress  history  is  also  reported  at  the  periphery  of  the  coneentrator  from  the  three  simulation  run 
solutions  (Fig.  14).  Initial  tensioning  of  the  polymeric  reflector  is  pronounced  at  each  projected  area  radius  of  all 
three  loading  scenarios.  These  stress  states  should  be  regarded  as  the  extreme  for  aperture  engagement  of  the 
respective  predictions.  The  subsequent  increase  in  distortion  energy  stress  manifests  upon  retensioning  due  to  a 
subsequent  excursion  cycle.  (Note:  the  simulation  run  solution  enforcing  the  step  wavefonn  body  flux  loading 
scenario  defined  by  Eq.  10  achieved  a  stress  state  relatively  higher  than  initial  tensioning,  with  respect  to  the  other 
two  scenarios,  as  the  velocity  for  this  loading  regime  was  opposite  sense  upon  concentrator  engagement,  returning 
from  an  excursion.  Only  the  initial  tensioning  stress  at  the  periphery  of  the  figure  acquisition  simulation  mn 
solution  enforcing  the  step  waveform  body  flux  loading  scenario  defined  by  Eq.  9,  valued  at  lO.lxlO’  Pa,  exceeded 
the  proportional  limit  of  Kapton  100  HN®,  valued  at  6.90x10^  Pa.)  Reasoning  for  the  conservatism  of  the  reflector 
stress  histories  due  to  model  reductions  is  supported  in  Ref  24. 


time,  s 

Fig.  14  MTR  subroutine  full  solution  of  concentrator  periphery  distortion  energy  stress  history  responding 

to  Eqs.  9-11  waveform  body  flux  loading  scenarios 


Reaction  forces  are  reported  as  the  sum  of  the  y-components  at  the  four  spine  vertex  nodes  (Fig.  15).  Nodes  of 
continuum  elements  do  not  have  rotational  DOFs;  reaction  moments  were  not  found.  (Note:  one  would  expect  the 
symmetric  x-component  forces  to  be  countered  by  the  opposing  spine  and  concentrator  sector  pairs  and  the 
asymmetric  and  net  y-component  forces  to  be  communicated  to  the  vehicle  bus.  Reaction  forces  at  the  model 
transition  are  truncated.  Maximum  x-  and  y-reaction  forces  for  the  three  loading  scenarios  were  approximately 
2,500  N  and  300  N,  respectively.) 


Fig.  15  MTR  subroutine  full  solution  of  spine  vertex  y-reaction  force  history  responding  to  Eqs.  9-11 

waveform  body  flux  loading  scenarios 


Temperature  history  of  the  reflector  periphery  can  be  plotted  throughout  the  entire  simulation  run  (Fig.  16). 
Thermal  authority  over  the  concentrator  was  found  locally  limited  to  the  periphery;  the  polyimide  conductivity  is 
from  one  and  two  orders  of  magnitude  less  than  the  conductivity  of  the  SMA  between  R  and  B2  crystalline  phases. 


Fig.  16  MTR  subroutine  full  solution  of  concentrator  periphery  temperature  history  responding  to  Eqs.  9- 

11  waveform  body  flux  loading  scenarios 


IV.  Conclusion 


A.  Deployment  Sequence 

The  deployment  sequence  modeling  efforts  were  insightful  for  both  forecasting  the  figure  acquisition  of  a 
gossamer,  polymeric  reflector  harnessing  the  SME  and  the  performance  capabilities  of  the  MTR  subroutine.  The 
three  waveform  body  flux  loading  scenarios  activated  R  to  B2  phase  transformations  of  the  SMA  spine  within 
approximately  177-374  s.  Distortion  energy  stress  and  temperature  histories  documented  through  these  simulation 


runs  are  indicative  of  the  mechanical  and  thermal  demands  on  the  two  continuums.  Features  of  these  load  trends  are 
relevant  to  other  SMA  actuation  schemes. 

Distortion  energy  stress  magnitudes  at  the  periphery  regions  of  the  spine  followed  the  phase  transformation 
kinetics;  greater  stresses  were  reported  midway  through  the  B2  crystalline  structure  recovery  corresponding  to  the 
more  rapid  transformation  rates  at  these  unstable  minima  and  consequently  greater  inertial  loads  communicated 
through  the  SMA  continuum.  The  mechanical  loads  were  significant  upon  engaging  the  aperture.  (Note:  these 
stress  spikes  are  within  Nitinol’s  austenite  proportional  limit  range. These  two  load  trends  are  reflected  in  the 
reaction  force  histories  of  the  simulation  run  solutions. 

The  distortion  energy  stress  magnitudes  at  the  periphery  region  of  the  reflector  reported  upon  initial  tensioning 
should  be  regarded  as  the  extreme  for  aperture  engagement.  The  subsequent  increase  in  distortion  energy  stress 
manifests  upon  retensioning  due  to  a  subsequent  excursion  cycle.  (Note:  the  simulation  run  solution  enforcing  the 
half-power  step  waveform  body  flux  loading  scenario  achieved  a  stress  state  relatively  higher  than  initial  tensioning, 
with  respect  to  the  other  two  scenarios,  as  the  velocity  for  this  loading  regime  was  opposite  sense  upon  concentrator 
engagement,  returning  from  an  excursion.  Only  the  initial  tensioning  stress  at  the  periphery  of  the  figure  acquisition 
simulation  run  solution  enforcing  the  full-power  step  waveform  body  flux  loading  scenario  exceeded  the 
proportional  limit  of  Kapton  100  FIN®.) 

These  mechanical  load  trends  decreased  from  the  full-power  step  to  the  step  waveform  body  flux  loading 
scenario  at  half-power  and  to  the  symmetrical  linear  waveform  body  flux  loading  scenario,  effectively  motivating 
the  phase  transformation  at  a  pace  tantamount  to  the  full-power  step.  Thermal  authority  over  the  concentrator  was 
found  locally  limited  to  the  periphery;  the  polyimide  conductivity  is  from  one  to  two  orders  of  magnitude  less  than 
the  conductivity  of  the  SMA  from  R  to  B2  crystalline  phases,  respectively.  All  three  body  flux  loading  scenarios 
imparted  an  equitable  amount  of  energy  to  the  system  and  consequently  the  final  temperature  realized  in  the  aperture 
near  the  interface  was  not  unique.  Again,  the  reflector  stress  histories  should  be  regarded  as  conservative  due  to 
model  reductions. 

Aspects  of  the  deployment  sequence  warranting  further  study  include  determining  the  effects  of  radiation  load 
paths  on  the  character  of  the  response  and  investigating  an  unreduced  model  to  perceive  true  three-dimensional 
phenomena  including  circumferential  stress  states  and  asymmetric  loading  of  the  reflector.  To  facilitate  figure 
maintenance,  detailed  models  could  incorporate  interfaces  between  the  spine  and  aperture  at  multiple  projected  area 
radii  and  account  for  imaging  and  electrostatic  hardware  affixed  to  the  spine. 

Additional  R  phase  configurations  should  be  considered  to  support  optimal  packaging  strategies  and  exact  B2 
phase  geometries  should  be  determined,  dictated  by  concentrator  position  and  reaction  force  demands.  Conceptual, 
collapsible  architectures  for  membrane  reflectors  of  merit  have  been  proposed.^^’  The  seminal  work  of  Mikulas^"' 
and  Fledgepeth^^  reasoned  large  space  vehicles  are  more  susceptible  to  on-orbit  loads  and  these  factors  drive  the 
basic  parameters  for  deployed  designs. 

In  addition  to  the  required  quasi-static  boundaries  afforded  by  the  spines,  concentrator  demands  include  stiffness 
to  reduce  the  wavefront  error  to  be  actively  managed  resulting  from  time  varying  disturbances.  Upon  reflector 
engagement,  the  integrated  model’s  first  excursion  cycle  frequency  excited  by  all  three  loading  regimes  was  non¬ 
unique,  regardless  of  the  temporal  character  preceding  the  aperture  tensioning.  This  tendency  suggests  a  coupling 
between  the  structural  dynamic  designs  for  deployment  and  in-service,  deployed  architectures.  Dynamic  stability 
scaling  laws  based  on  structural  modal  frequency  and  damping  ratio  relating  aperture  dimensions  to  structural  depth 
and  material  specific  modulus  have  been  presented  by  Peterson  and  Flinkle.^®  Structural  deformation  sensitivity  was 
found  to  scale  inversely  with  the  square  of  structural  depth  and  inversely  with  the  material  specific  modulus. 

Parallel  research  complimentary  to  the  spine  and  concentrator  model  has  been  performed  concerning  the 
management  of  in-service  aberration  phenomenon.  The  scope  of  this  study  consists  of  characterizing  the  non-linear 
dynamic  behavior  of  membrane  reflectors  to  visible-optics  quality  under  realistic  support  and  loading  scenarios. 
The  objective  of  this  work  was  to  investigate  the  non-linear  characteristics  of  the  dynamic  behavior  of  membrane 
optics  as  the  basis  for  future  work  in  system  identification  and  figure  control.^’  Further  considerations  of  model 
error  limits,  constitutive  relation  derivations,  and  numerical  implementation  have  been  combined  by  deBlonk^*  into 
a  guideline  for  modeling  the  optical  level  structural  behavior  of  curved  membrane  mirrors. 

Continued  development  of  the  MTR  subroutine  to  support  two-way  SME  behavior  would  be  of  interest. 
Scalability  of  the  parameter  window  mapped  for  one-way  behavior  from  replicating  experiments,  with  SMA  test 
articles  an  order  of  magnitude  less  in  size  and  undergoing  relatively  less  detwinned  strain  recovery,  deserves  some 
scrutiny.  It  is  evident  from  stills  of  the  spine  FEMs,  logic  implementing  compatibility  theory  to  ensure  conservation 
of  mass  is  required  for  greater  fidelity  of  SME  predictions  undergoing  an  R  to  B2  phase  global  geometry  return, 
with  relatively  large  SME  attendant  translations  and  rotations. 


B.  MTR  Subroutine 

The  development  and  refinement  of  an  SMA  modeling  approaeh,  the  MTR  subroutine,  has  been  presented.  In 
eoneert  with  the  eommereial  finite  element  solver,  ABAQUS  Explieit,  this  applieation  supports  multi-dimensional, 
temperature-displaeement  transient  predietions  of  the  SME  exhibited  by  SMAs  through  implementation  of  a 
phenomenologieal  eonstitutive  law.  The  subroutine  eapability  is  limited  to  one-way  phase  transformation,  from 
martensite  to  austenite,  attendant  detwinned  strain  reeovery;  two-way,  superelastieity,  and  eonventional  plastie 
deformation  behavior  is  not  supported. 

Efforts  to  tune  the  approaeh  against  SME  experiments  foeused  on  mapping  parameters  governing  the  SIM 
phenomenon  and  struetural  damping.  Due  to  eomputational  eapaeity  limitations,  the  MTR  subroutine  was  retooled 
to  interfaee  with  the  ABAQUS  Standard  implieit  direet  integration  solver.  Numerieal  eomplieations  were 
eneountered  with  attempts  to  replieate  seenarios  involving  eonstant  natural  boundary  eonditions. 

Although  further  advanees  on  both  governing  parameter  seores  is  reeommended,  model  refinement  was 
eonstruetive  for  identifying  the  undesirable  numerieal  artifaets  whieh  arise  when  modeling  eonstant  natural 
boundary  eonditions.  For  seenarios  not  involving  eonstant  loads,  a  published,  experimentally  measured  non- 
dimensional  value  of  this  Nitinol  material  property  was  used.^^  Also,  insight  was  gained  regarding  insensitivity  to 
the  struetural  damping  governing  parameter  at  the  experimental  phase  transformation  rates;  a  revamped 
experimental  design  would  be  required  to  map  this  parameter.  In  addition  to  mapping  the  SIM  parameter  window, 
refinement  of  the  phase  transformation  kineties  model  with  the  addition  of  a  series  of  eorreetion  terms  to  eapture  the 
features  of  the  SME  distribution  of  a  polyerystalline  SMA  aetuator  was  of  value. 

The  greatest  asset  of  model  refinement  is  the  eonfidenee  established  knowing  doeumented  numerieal  instabilities 
and  divergenee  trends  of  the  MTR  subroutine  modeling  approaeh  are  not  assoeiated  with  simulation  mns  involving 
inertial  loads.  This  eredenee  staged  the  foreeasting  of  the  transient  response  of  a  spaeebome  membrane  optie  SMA 
deployment  aetuator  spine  system,  the  headlining  objeetive  of  this  thesis  researeh,  and  also  establishes  eonfidenee  in 
the  amenability  of  the  approaeh  to  sehemes  implementing  the  SME  to  reeonfigure  and  tune  monolithie  struetures  or 
aetuate  meehanieal  systems  not  involving  eonstant  natural  boundary  eonditions. 

Aspeets  of  the  MTR  subroutine  warranting  further  study  inelude  traeking  and  quantifying  the  systemie 
influenees  of  the  boundary  eondition  steps’  energy  eontributions  and  timing  of  the  simulation  run  step.  The  effeet 
order  on  the  simulation  run  solution  of  assuming  the  SIM  angle  as  a  eonstant  value  through  the  phase  transformation 
is  worth  eonsideration  to  perhaps  enable  greater  fidelity.  Pertinent  issues  arising  from  the  deployment  sequenee 
entailed  the  sealability  of  the  parameter  window  mapped  from  replieating  experiments,  with  SMA  test  artieles  size 
and  relative  detwinned  strain  reeovery,  and  the  inelusion  of  logie  implementing  eompatibility  theory  to  ensure 
eonservation  of  mass  of  shape  ehanging  eontinuums. 

Continued  development  of  the  MTR  subroutine  to  support  two-way  SME  behavior  is  a  manifold  issue  for  taekle. 
Revising  the  approaeh  to  support  this  phenomenon,  would  render  the  subroutine  eapable  of  throughput  paekaging 
and  deployment  design.  One  eould  envision  a  reversal  of  the  eonstituent  hypothesis  to  permit  and  aeeount 
detwinning  strain  aeeumulation  at  the  threshold  stress  and  until  the  maximum  reeoverable  strain  limit. 

Computational  eapaeity  was  a  reeurring  theme  throughout  this  investigation.  Resouree  limitations  preeluded 
interfaeing  with  the  ABAQUS  Explieit  solver,  as  the  MTR  subroutine  was  designed,  and  a  dynamie,  fully-eoupled 
thermal-stress  analysis  ineluding  inertia  effeets  were  forgone.  Reeognizably,  models  foreeasted  were  spatially  and 
temporally  eoarsely  meshed  to  minimize  eomputational  expenses.  Enhaneements  in  eomputational  eapaeity  will 
serve  to  expand  the  feasibility  of  this  SMA  modeling  approaeh. 

Just  as  spaee  vehieles  first  ineorporated  the  SME,  deployment  and  in-serviee  ehallenges  of  gossamer  spaeeeraft 
sueh  as  the  transient  response  of  a  spaeebome  membrane  optie  SMA  deployment  aetuator  spine  system  addressed 
through  this  thesis  researeh,  will  further  showease  the  leverage  SMA  teehnology  extends  to  straetural 
reeonfiguration  and  tuning  to  foster  and  mature  sueeessive  large,  launeh-paekaged  spaee  vehiele  programs. 
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